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1. Introduction 

The effect of disorder in high-T c cuprates, either intentionally introduced or present as a 
result of the doping process, has always been crucial to their practical applications. For 
several classes of cuprates, Scanning Tunneling Microscopy(STM) has been particularly 
useful to clarify how disorder affects superconductivity [U El El IU El E] . In the frequently 
studied Bi 2 Sr 2 CaCu 2 08+ ( 5(BSCCO) system, where disorder is induced primarily by 
interstitial oxygens which dope the system, STM data show a very pronounced 
inhomogeneity of superconducting gap over the sample surface. In particular, the 
gap magnitude is significantly enhanced in the vicinity of oxygen defects pE]. A 
phenomenological model was proposed by Nunner et al [7J based on the assumption 
that the pairing interaction was locally enhanced in the vicinity of oxygen defects, 
which successfully captured the feature of locally enhanced superconducting gap, and 
the reduction of coherence peak in the large gap regions, as well as other statistical 
correlations. These authors pointed out that such an enhancement could be understood 
from a strong coupling perspective if the local superexchange J(r) was increased by 
the presence of a nearby dopant. Maska et al [8] then showed that within a single- 
band Hubbard model, local variations in site energies as might be due to a dopant 
perturbation would always lead to an enhanced J within second-order perturbation 
theory, but similar calculations within a three-band Hubbard model [9] indicated that 
the sign of J modulation was not universal. Johnston et al [10] found a suppression 
of the superexchange in Cu-0 cluster calculations with point charge estimates of 
Madelung energies, but pointed out that enhancement could occur if the electron- 
phonon coupling to c-axis oxygen modes were accounted for. From a weak-coupling 
perspective, Foyevtsova et al [llj considered the change in local electronic structure 
due to an interstitial O in first principles calculations of BSCCO, and used the results 
to calculate the spatially averaged spin fluctuation interaction in undoped and doped 
system; the enhanced pairing in the doped system was similarly attributed to local 
enhancement around O interstitials. A somewhat different approach has been pursued 
by other groups [T2J, [13] , who have discussed the large polarizability of the O interstitial 
as possibly modulating the overall polarizability of the Cu-0 complex and hence the 
pairing. 

Although the inhomogeneous superconducting gap and several important 
correlations among STM observables are nicely reproduced by the locally enhanced 
pairing interaction picture at optimal doping, a puzzle remained: if each oxygen defect 
locally enhances superconducting gap, why does the average gap decrease with increasing 
oxygen concentration? This feature has been quantified in [5], where empirically not 
only is the average gap found to be decreasing with oxygen concentration, but the 
spatial distribution of superconducting gap is observed to be narrower at large doping. 
This apparent discrepancy indicates that additional ingredients are necessary in order 
to understand the doping dependence of the inhomogeneous superconducting gap. A 
natural assumption is that at least one of these additional ingredients must be global 
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electronic correlations, which affects the average superconducting gap over a rather large 
scale, but still allows the local gap to be enhanced around each oxygen defect. 

In this work we employ slave boson mean field theory to study the inhomogeneous 
pairing gap on the surface of BSCCO, and show that besides the locally enhanced 
pairing interaction, the holes donated by oxygen defects have a dramatic effect on the 
pairing gap. We find that at large oxygen concentration, the hole doping level increases 
accordingly, which reduces the density of states(DOS) of spinons and causes the average 
gap to drop. The advantage of slave boson mean field theory is that it treats disorder and 
strong correlations on equal footing, since the inclusion of strong correlations is necessary 
to explain the reduction of DOS. A theory of this kind has been applied to study 
YBa 2 Cu 3 6+x .(YBCO) in the presence of a point-like repulsive impurity [H], where the 
repulsion of holes and the reduction of the fermionic bandwidth in the vicinity of the 
impurity influence charge and spin response significantly, yielding a Nuclear Magnetic 
Resonance(NMR) spectrum consistent with experiments. In essence, oxygen defects 
in BSCCO serve as impurities of the opposite kind: their Coulomb potential attracts 
holes which, in the slave boson description, promotes the enhancement of the fermionic 
bandwidth. We will demonstrate that the enhancement of fermionic bandwidth serves 
as another mechanism to reduce the coherence peak in the large gap regions. 

A similar method [15] has been proposed previously to study the effect of oxygen 
defects in BSCCO, where screening of the ionic Coulomb potential is studied in detail, 
and is argued to be the origin of electronic imhomogeneity. However, without a locally 
enhanced pairing interaction, the superconducting gap is found to be anticorrelated 
with oxygen defects, in contrast to the observation of STM. It was also proposed that 
postulating two types of oxygen dopants, one playing the role of Coulomb defect and the 
other purely as controlling the doping, can resolve the O-gap anticorrelation issue [16J. 
However, based on various microscopic calculations [121 1131 El El [10] that support 
the possible enhancement of the local pairing interaction by oxygen defects, we believe 
that it is important to include this effect in a phenomenological way, and treat each 
oxygen defect as both a Coulomb center and a source of charge. We will show that the 
positive correlation between pairing gap and oxygen defects is recovered by properly 
accounting for various energy and length scales in the problem. It is also proposed that 
slave boson treatment naturally enhances the pairing interaction around oxygen defects, 
without any additional microscopic mechanisms such as spatially varying atomic level 
[T7] . The argument is based on the renormalization of hopping integral and its 
relation with the superexchange interaction Jy = 4tfj/U. However, this mechanism 
gives a very small enhancement of Jy, since we found that the hopping integral is only 
renormalized by about 20% near oxygen defects, which cannot account for the large gap 
inhomogeneity observed in STM. The inclusion of a phenomenological, locally enhanced 
pairing interaction in the framework is therefore still necessary in the development of 
the theory. 

The structure of the paper is as follows: in Sec. II we revisit the enhanced pairing 
interaction picture in [7], and study the change of local bandwidth due to strong 
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correlations phenomenologically. This section essentially illustrates all the physical 
effects that occur in the vicinity of each oxygen defect. In Sec. Ill we introduce the slave 
boson mean field treatment for homogeneous cuprates, and demonstrate the well-known 
reduction of average gap due to hole doping. In Sec. IV, oxygen defects are introduced 
into the slave boson treatment to study disorder and correlations on an equal footing. 
Sec. V concludes with our results. 

2. phenomenological model of a single oxygen defect 

In this section, we follow the strategy of [7] to study a single oxygen defect within a 
Hartree-Fock-Gor'kov mean field approach. The motivation is to clarify the influence of 
oxygen defects on the local density of states(LDOS), in the simplest non-self-consistent 
manner. In comparison with the fully self-consistent treatment in Sec. IV, only the 
fermionic degrees of freedom are considered in this section, since they are essential for 
local pairing, and the effect of strong correlations is discussed phenomenologically. We 
consider the following impurity Hamiltonian 

ijcr icr 

+ £ (Atf + dA ij )c^c] i + h.c. (1) 

For simplicity and consistency with [TJ [18] , in this section we assume that the oxygen 
defect is located on top of a random Cu site (this is not quite correct [E]), and we 
choose t = 1 = 400meV for nearest-neighbor and t' = —0.3 for the next-nearest- 
neighbor hopping in this section. The homogeneous pairing amplitude is Ay = 
Ao (5i =r 5j =r ± x — 5i =r 5j =r ± y ) . In [7], it is also shown that if only the electronic (fermionic) 
degrees of freedom are accounted for, the usual association of an on-site Coulomb 
potential V* to each oxygen defect gives a form of the LDOS, including a highly 
inhomogeneous low-energy LDOS, that is inconsistent with STM experiments. Here 
we focus on the effect of hopping impurities dtij and the off-diagonal (pairing channel) 
impurities dAij, and assume that they only differ from the bulk value on the four bonds 
connected to impurity site r, dtij — dtji — dt (5i =r 5j =r ± Xtr ± y ), and dAy = dAji = 

dA (Si =r Sj =r ± x 3i=r3j=r±yJ • 

For comparison with [7J [18] , we first perform simulations for a single off-diagonal 
impurity. As shown in figure Ufa), one clearly sees a enlarged gap and reduced coherence 
peak for the case of local pairing enhancement dA/ ' A$ = 1. The off-diagonal impurity 
case obviously captures the enhancement of the gap by oxygen defects. We then perform 
simulations for the hopping impurity, for both hopping enhancement and reduction 
cases. The motivation comes from the slave boson treatment of BSCCO surface, where 
the spinon hopping amplitude is enhanced near oxygen defects due to accumulation of 
holons, as will be addressed in Sec. IV. Here we treat this effect at the phenomenological 
level by using a hopping amplitude locally modulated by roughly dt/t ~ 10% as found by 
the microscopic theory discussed in Sec. IV. Figure [TJ(b) shows the simulation for both 
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Figure 1. LDOS on the impurity site of the phenomcnological model described 
by equation (|T|), with different impurity models: (a)Pairing enhancement /reduction 
dA/Ao = ±1 on four bonds connected to the impurity site, (b)hopping 
enhancement /reduction dt/t — ±0.1 on four bonds, and (c)both pairing and hopping 
enhancement / reduction. 



hopping enhancement and reduction on four bonds connected to the oxygen defect. 
The LDOS on the impurity site clearly indicates that hopping enhancement reduces 
the coherence peak and vice versa, while the gap size remains unchanged. The LDOS 
calculated over a wider energy range, as shown in figure [T^b), demonstrates that the 
reduced coherence peak is due to additional bound states formed on the top and the 
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bottom of the band. Since the overall spectral weight has to be conserved, the formation 
of bound state reduces the spectral weight away from the band edge, particularly at the 
coherence peak. 

The simulation that contains both pairing and hopping enhancement is shown in 
figure [0(c), where one sees that the superconducting gap is enlarged as expected, and the 
coherence peak is further reduced compared to the case that contains only an off-diagonal 
impurity. On the other hand, if the impurity locally reduces both the pairing and 
hopping interactions, then the gap is reduced and the coherence peak height increases. 
This simulation suggests that the change in the local hopping amplitude, whose origin 
can be traced back to the holon accumulation due to strong correlations, plays an 
important role in the shape of STM spectrum, particularly the weight of coherence peak. 
The calculation that associates change of hopping amplitude to correlation effects will 
be demonstrated in Sec. IV. 

3. slave boson mean field theory in a homogeneous cuprate 

Before we introduce the disorder model for BSCCO surface in the next section, we 
present here the slave boson mean field model of a homogeneous cuprate. One can 
regard the phenomenological impurity model in Sec. II as a demonstration of how 
oxygen defects affect STM spectrum locally, while the mean field theory in this section 
demonstrates how oxygen defects change the average gap value globally. In the slave 
boson mean field treatment at finite doping, spin and charge degrees of freedom can 
be separated but mutually renormalize each other. The pairing gap is described by the 
pairing of spinons, which relies on the LDOS of spinon band that in turn is sensitive to 
hole doping. 

We begin with the t — if — J model treated with the slave boson method 



We choose t = 400meV= 1 and next-nearest-neighbor hopping t' = —0.4 , a slightly 
different value from Sec. II such that the van Hove singularity always locates close to 
but lower than the pairing gap. The exchange coupling J = 72meV= 0.18 is smaller 
than the well accepted value J f=s 130meV, hence one expects a smaller gap in the 
homogeneous limit. However, we found that the local enhancement of J due to oxygen 
defects restores the pairing strength, and results in a pairing gap comparable to the 
STM data, as will be demonstrated in Sec. IV. 

In equation fl2]), the electron operator has been decomposed into the spin(spinon) 
and charge(holon) component, c lcr = f% a b\. n{ a = f\ a fi a and n\ = b\b{ represent 
the spinon and holon density, respectively. In the homogeneous case, the nondouble 
occupancy constraint J2a flaf™ + b\h = 1 can be satisfied by adjusting spinon and 





(2) 
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holon chemical potential independently, so the Lagrangian multiplier A, is practically 
irrelevant. We apply the following mean field decoupling to the Heisenberg term 

Si ■ S 3 = SIS* - - } I ),,!)„ + W nlnf_ a . (3) 



where 



Si — ^fiaVapfip , 



(4) 



The Hamiltonian is then 

n = £ -*u hb]flf 3 , + £ j (s-s* - \d\ 3 d % 

ija (ij) 
ia i icr 

(5) 

Notice that the spinon density-density term nfn|/4 is cancelled. The resulting 



Hamiltonian is slightly different from other slave boson decoupling schemes [20, 21J. 
However, we do emphasize the importance of keeping the SfSj term. A mean field 
decoupling of this kind shows that S?Sj accounts for the magnetic correlations that 
cause impurity induced magnetization in the normal state [H]. Evidence from neutron 
scattering indicates an induced magnetization in the center of a vortex [22J, which 
suggests that magnetic correlations also exist in the superconducting state. We therefore 
choose to adopt the mean field decoupling in equation ([3]), as it is crucial to other 
calculations involving magnetic field. 

Renormalization of the holon bandwidth is associated with the following operator 

Kija — ficrfja i (6) 

while the renormalization of the spinon bandwidth is simply associated with hi at low 
temperature. In the mean field treatment, operators in equations (jl]) and are replaced 
by their ensemble average 

(Sf) = tr* , 

(4) = d *i . 

(&I> = % , 

(blbf) = Qij . (7) 

The nonzero value of (b\) in equation ([7]) implies Bose-Einstein condensation(BEC) in 
the holon sector. At temperature relevant to the STM measurement that we primarily 
compare with [5], roughly T$tm < 50K, we find that the holon ground state takes more 
than 99% of the hole population at any doping larger than 5%. Thus one can safely 
assume that bosonic order parameters are entirely determined by ground state wave 
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functions, and we expect (bjbj) = (b\)(bj) = (see below). The usual mean field 
procedure is applied to the above formalism [21], which yields the following Lagrangian 
for spinons and for bosonic excitations 



+ £ Ai ( £ fUia + - 1) + E ^ 

i \ <y ! (ij) 

+ E ~tijQijfiafj<r + E ( _ ^ E ) ^7 

ij'cr ij i \ / 

+ E-^A. + ^.c. (8) 

In the homogeneous d— wave superconducting state, the following ansatz is made for the 
order parameters 

_ * _ o 

X{ij)tr — X{ij)a ~ X i 

_ * _ 0/ 

X{(ij))a - X({ij))a - X , 

Q(ij)a = Q*ij)a — Q° ) 
Q({ij))a = Q*{ij))a = Q°' > 

d - = d* = { d ° f ° r J = i ± & 
13 13 \ -d for j = i ± y 

mi = 0, (9) 

where (ij) and ((ij)) denote nearest- and next-nearest-neighbor bonds, respectively. 
The Lagrangian of holon and spinon are then separated. The spinon part has dispersion 
and d— wave gap given by 

£k = —2tQ° (cos k x + cos k y ) — At'Q ' cos k x cos k y — fj,f , 

Ak = Jd (cos k x — cos k y ) , 



k \ 



(10) 

For the holon part, The Hamiltonian is 

k 

Wk = — 4tx° (cos + cos k y ) — 8t'x°' cos k x cos k y — ^ . 

(11) 

The order parameters are solved self-consistently by means of Bogoliubov-de 
Gennes(BDG) equation with the inclusion of BEC of holons. For the spinon sector, 
we apply the following spin-generalized Bogoliubov transformation 

fit = E M «,iT7nT + <it^4 
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fil = Wn,47r4 + <,4Tnt ■ ( 12 ) 
n 

After diagonalization, the spinon order parameters are calculated by wave functions 
u n ,i<j and f n ,iCT) and Fermi distribution f(E na ) 

(flfja) = { U *n,ia U n,jJ(Ena) 
n>0 

+ K,i* V n,jo(l ~ f(Kw))} , 

(fafi+8l) = { U n,i\ V *n,i+5i{ 1 ~ f{ E n\)) 

n>0 

+v:^u nA+Si (l-f(-E ni ))} . (13) 

Bosonic order parameters under BEC are calculated by assuming all holons are 
condensed in the ground state 

(h) = v/iVoao,, , (14) 

where a 0i { is the ground state wave function of equation ffTTj) . and N = xL 2 is 
total number of holes on an L x L lattice. The d— wave superconducting state is 
characterized by the region in the phase diagram where 4 ^ and T < Tbec- 111 
reality, BEC in slave boson formalism is a rather complicated issue. As pointed out 
in [21], if one naively adopts the BEC formalism for 2D interacting bose systems, the 
condensation temperature is as high as T BEC ~ O(1000)K near optimal doping, which 
is certainly unrealistic. Such a high condensation temperature is properly reproduced 
if one generalize the BEC treatment, equations ([7]) and ( IT4|) to finite temperature. To 
reduce the condensation temperature to a realistic scale, one needs to include the gauge 
degrees of freedom and their scattering with holons [21] . Nevertheless, we are primarily 
interested in low temperature STM results deep inside the condensate, T <C Tbec- 
Therefore the precise value of Tbec does not affect the results in the present work, and 
our assumption that holons occupy only the ground state, equation ( !T4|) . is adequate. 
For the calculation of gap inhomogeneity, we choose T = 50K such that it sits reasonably 
below T c in most of the doping regime, while safely above the temperature limit due to 
finite cluster size. STM data at this temperature is also available to compare with. 

Figure |2] shows the pairing gap and DOS as a function of hole doping. Notice that 
since holons remain condensed, STM effectively measures the spinon spectrum. The 
implication of figure [2] is straightforward: as doping increases, the spinon bandwidth 
enlarges, causing reduction of DOS at Fermi surface, hence a smaller pairing gap. 
The enlargement of the spinon band is due to a larger Q Q . This is because at low 
temperatures, Q m n b = x, so the spinon bandwidth is roughly proportional to doping. 
Our simulation simply indicates that, although the pairing interaction is enhanced 
locally around each oxygen defect, the fact that oxygen defect donates holes into Cu02 
plane can change the spinon DOS dramatically. This serves as another mechanism to 
determine the pairing gap in a larger length scale, as long as holes are practically not 
localized. The precise simulation that contains both local and global effects is presented 
in the next section. 
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Notice that in order to compare our simulation with real cuprates, the hole doping 
level has to be rescaled. It is known that slave boson mean field theory of this kind have 
problems getting both temperature and doping scale correctly [23j [21] . We choose our 
parameters such that the simulation yields a pairing gap close to the value indicated 
by STM(A ~ 20meV at optimal doping). The doping level is then rescaled empirically 
by identifing the doping where superconductivity vanishes x ~ 0.6, with that of a real 
cuprate p fa 0.3, therefore 

x fa 2p. (15) 

Throughout the article we denote numerical doping scale as x, and the true doping 
scale in cuprates as p. One should rather regard equation (115j) as an empirical 
formula, since the precise correspondence between x and p depends on the mean field 
decoupling equation ((3]). Nevertheless, the pairing gap and critical temperature at 
optimal doping(T c fa 100K at p = 0.15 = x/2) is very close to the true value in BSCCO. 
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4. slave boson mean field treatment of BSCCO surface 

In this section, we incorporate oxygen defects into slave boson mean field theory studied 
in Sec. Ill to simulate the surface of BSCCO. We adopt the phenomenological treatment 
of [7] (see also [17]) that assumes oxygen defects locally enhance pairing strength and 
attract holes to their vicinity. In contrast to the weak coupling theory in Sec. II, 
enhancement of the pairing interaction is described by locally varying in the t — t' — J 
model, 

Jij = J + SJ(h i + h j ) /2 , (16) 

and the Coulomb interaction yields an additional term in the Hamiltonian described by 
locally varying potential Vi 

^S eXp( - r " /A> , (17) 

s r is 

where s is the projected position of the oxygen defects into the plane. The precise value 
of S J, Vq, and A will be addressed later in the discussion of different energy and length 
scales. With all these effects together, the Hamiltonian in the slave boson language 
reads 

H f = J]] ~tijQijfi<jfja + ~tijQijfjafi<J 
ija ija 

+ £ j i3 ,s; mj + £ - J f dl 3 D« + £ --fd^ 

(ij) (ij) (ij) 



H b = -tijXijablbj + Yl -tijXijab]bi 



ija ija 
b 



n 



Notice that Vi is present only in the bosonic sector, as it only affects charge degrees 
of freedom in the problem. In the disorder case, Aj serves as a local correction to the 
chemical potential in order to satisfy the nondouble occupancy constraint. To compare 
with the density inhomogeneity in [23], we assume that oxygens are located on top of 
four adjacent Cu sites, although this is slightly different from the most energetically 
favorable position [19]. We assume the positions of oxygen defects are truly random, 
unlike the nonadjacent condition employed in [25] • This is because if the oxygens are 
assumed to be nonadjacent to each other, then in a 2D square lattice their concentration 
can at most be 25%, whereas a higher oxygen concentration is necessary to achieve higher 
doping level in the present work. 

Figure [3] shows the probability distribution P(A) and real space map of pairing gap. 
In figure [3][a), one sees a fairly broad distribution for the doping regime 0.1 < x < 0.5, 
whose mean shifts to smaller value as doping increases. This is a clear indication that the 
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Figure 3. (a) Probability distribution P(A) of the pairing gap at different doping, 
in a 30 x 30 cluster averaged over 40 impurity configurations. (b)The average gap 
and (c)the standard deviation calculated from -P(A). (d)~(g) Real space gap map 
for a specific impurity configuration in different doping levels. Oxygen positions are 
indicated by white dots. The temperature is fixed at T = 50K in this figure. 

average pairing gap decreases with increasing doping, consistent with the observation in 
STM. Comparing figure E^b), where the average pairing gap (A) is quantified for each 
probability distribution, with the homegeneous gap shown in figure El it is evident that 
this (A) versus x behavior is inherent from the doping dependence in the homogeneous 
case. Since the lineshape of -P(A) is rather peculiar, we calculate its linewidth by means 
of standard deviation, as shown in figure [3^c) . In comparison with the STM data [5], 
where the linewidth narrows significantly, about 60% from underdoped to overdoped 
regime, our simulation shows a narrowing only about 20%. We suspect that to correctly 
capture this 60% reduction of -P(A) linewidth, one needs a fine mechanism to control the 
overlap between superconducting islands as doping changes, for instance a fine tuning of 
A, or using a more sophisticated impurity model other than equation (JTBl) and equation 
( I17p . Nevertheless, our calculation demonstrates that within a simplestic impurity model 
with only one length scale A associated with impurity, we can qualitatively capture the 
pairing gap inhomogeneity induced by oxygens. Notice that the linewidth of -P(A) 
increases again in highly overdoped case x = 0.5. This is because in the absence of 
oxygen disorder, the homogeneous gap is zero at this doping. Randomly distributed 
oxygens locally enhance the gap, which causes the distribution to range broadly from 
zero to A w 0.08. 

From the real space gap map, figure [3]^d) to (g), the inhomogeneity induced by 
oxygen defects (white dots) is clearly visible. For each doping level, the pairing gap is 
locally enhanced near oxygen defects, consistent with the weak coupling approach [?]. 
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This is particularly evident at low doping level, as shown by red spots in figure El^d) for 
x = 0.1, where large gap region coincides with the oxygen positions. The oxygen-rich 
and oxygen-poor regions also give a double line structure to the gap distribution shown 
in figure [3]^a), particularly for low doping lines from x = 0.1 to 0.3. At higher doping, 
puddles of large gap region start to overlap, which eventually smears out the double line 
structure at around x = 0.4. 

We now discuss the robustness of the simulation against variation of the parameters, 
especially interaction strength 5J and Vo, and interaction range A. We find that in 
order to have an average pairing gap decreasing with doping, A has to be relatively 
small, roughly A < 1. In other words, the interaction has to be modulated on an 
atomic scale j7]. This can be understood simply because, had the pairing enhancement 
due to oxygen defects been longer range, it would have overcome the correlation-induced 
doping effect and enhanced the gap globally. Secondly, in order to produce a probability 
distribution that has comparable linewidth to STM data, the pairing enhancement must 
be roughly the same as the homogeneous superexchange interaction 5 J ~ J. This gives 
an effective coupling 1.34J on the nearest Cu-Cu bonds of an oxygen defect. The 
strength of Coulomb interaction Vo is also crucial to pairing gap. We find that if Vo is 
too large, the holon density is highly accumulated near oxygen defects, causing a very 
low spinon density and a smaller gap, even though pairing interaction is enhanced. As 
also found by Nunner et al., a large value of Vo can destroy the anticorrelation between 
the coherence peak width and position. This analysis sets up an upper limit for the 
Coulomb interaction Vo < 0.4. We choose Vo = 0.2, close to the value at optimal doping 
obtained by estimation of the dielectric constant |25j, although the value here is taken 
to be doping independent. 

A simulation of the LDOS and order parameters for a specific disorder configuration 
is shown in figure HJ For any doping level examined, the height of the coherence peak is 
found to decrease with gap size and vice versa, coinciding with [7] and the STM data. 
Comparing the Aj,- map in figure H^f) with the Q^ map in figure Wig), it is evident that 
regions with large gap also have large Qij. Similar to the homogeneous case, this is 
because at relatively low temperatures, bosons are fairly condensed, so ~ ^jn\n y 
Thus is large near oxygen defects because the holon density is rich. Since Qij 
enters the spinon Hamiltonian equation ( fl8l) as hopping amplitude, enhancement of 
Qij means the spinon hopping amplitude is locally enhanced around oxygen defects, 
as assumed phenomenologically in Section II. On the other hand, the LDOS near the 
oxygen defects has a very pronounced impurity bound state at band edge, which further 
reduces coherence peak. From the phenomenological model in Sec. II, it is clear that 
this bound state is a result of hopping enhancement, i.e. larger Q^ near oxygen defects. 
Therefore from this calculation of the LDOS, one can see the importance of strong 
correlations on local bandwidth and how it influences the shape of LDOS. 

Figure Hie) shows a very significant holon density inhomogeneity. This compares 
well with [25], where inhomogeneous charge density on BSCCO surface was studied by 
means of small hole pockets, For average doping x = (nf) = 0.33, the holon density varies 




Figure 4. (a)^(c)LDOS at three different average doing. In each figure, solid, dashed, 
and dotted lines represent typical LDOS at positions with large, medium, and small 
gap, respectively. Inserts show the LDOS near chemical potential. (d)~(g)Local 
order parameters at a specific impurity configuration at x = 0.33, in a cluster size 
7.6nmx7.6nm. 



from 0.3 < n\ < 0.4. However, in the weak coupling mean field approach of [?J, the hole 
density inhomogeneity is very moderate, roughly less than 5%. This can be understood 
as the following: in both [25] and the present work, charge degrees of freedom are carried 
by dilute holes, which screen Coulomb impurities poorly, while in [7] it is carried out 
by electrons, which screen Coulomb impurities very efficiently We anticipate that an 
approach like the present work is particularly accurate for underdoped to optimally 
doped cuprates, where hole density inhomogeneity is due to poor screening of dilute 
holes. On the overdoped side, where the system has presumably crossed over to a Fermi 
liquid, one should perhaps adopt weak coupling approaches like in [7] to accurately 
describe the screening by electrons. Notice that the density inhomogeneity in the present 
work is slightly smaller than in [25], due to coupling between charge and spin degrees of 
freedom. This is similar to a previous report where no clear Friedel density oscillation 
is seen near a point-like repulsive impurity [H]. This is simply because holons, which 
respond to the Coulomb interaction of disorder, do not have a characteristic length scale, 
so its density profile is rather smooth. Due to the nonodouble occupancy constraint, the 
amplitude of charge inhomogeneity is also suppressed, because holon density is subject 
to the spinon density that does not respond to Coulomb interaction directly. 




1 2 3 4 5 6 7 

R/a 



Figure 5. Various correlation functions calculated according to [25], at doping level 
x = 0.3 and T = 50K. 

Lastly, we address the issue of various correlation functions related to STM. We 
define the correlation functions according to [25], which is consistent with [3] but defined 
on a discrete lattice. The results for various correlators are shown in figure [5j The oxygen 
defects are positively correlated with pairing gap. The value of correlation extracts to 
~ 0.7 for onsite correlation R = 0, close to the value ~ 0.4 extracted from STM [4]. The 
holon-oxygen correlation, which is exactly the negative of spinon-oxygen correlation due 
to the no double occupancy constraint, is always positive, indicating the fact that an 
oxygen defect attracts holes to its vicinity. The length scale at which correlation drops 
to zero is roughly 2a, which is smaller than that reported in [25]. Notice that in [TJ, 
the potential length scale A = 0.5 yields a correlation length in good agreement with 
the STM results, whereas the same value of lambda gives a smaller correlation length in 
the current approach. We suspect that this is because spinon correlation functions are 
subject to holon correlation functions, and the no double occupancy constraint reduces 
the correlation length. Since A is treated phenomenologically, in principle it can be 
adjusted to yield a longer correlation length. However, in the present approach - where 
the oxygen concentration is rather high - our choice of a smaller lambda allows us to 
keep the pairing enhancement local. 

5. Conclusion 

We have presented a slave boson mean field approach to compare the possible physical 
effects induced locally by oxygen dopants near the surface of BSCCO, and show how 
their interplay can explain the pairing gap revealed by STM experiments. Firstly, the 
enhancement of the pairing gap and reduction of coherence peak near oxygen defects 
is attributed to locally enhanced pairing interaction via superexchange. Secondly, the 
average spinon density of states at the chemical potential is reduced at large doping, 
causing the average gap to decrease, although it remains locally enhanced around oxygen 
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defects. We also found that since holes are attracted to the vicinity of oxygen defects, 
an impurity bound state is formed at the edge of spinon band, which further causes the 
reduction of coherence peak, yielding a LDOS shape consistent with STM. 

The present work also gives important evidence regarding various length scales and 
interaction strengths in the cuprates. Since the pairing gap is locally enhanced around 
oxygen defects but has an average value that decreases with doping, we conclude that the 
length scale at which the pairing interaction is enhanced is very small A ~ 0.5a, although 
the enhancement 5J is about the same magnitude as the bare pairing interaction 
J, such that the it gives a significant gap inhomogeneity. On the other hand, the 
Coulomb interaction associated with the oxygen defects is relatively small compared 
to the hopping amplitude, such that the accumulation of holes around oxygen defects 
does not lead to a reduction of the pairing gap near oxygen defects, although it is large 
enough to yield a significant charge inhomogeneity in underdoped regime. 
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